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Abstract 

This paper is concerned with dynamic user equilibrium (DUE) with elastic travel demand 
(E-DUE). We present and prove a variational inequality (VI) formulation of E-DUE using 
measure-theoretic argument. Moreover, existence of the E-DUE is formally established 
with a version of Brouwer's fixed point theorem in a properly defined Hilbert space. The 
existence proof requires the effective delay operator to be continuo us, a regularity con dition 
also needed to ensure the existence of DUE with fixed demand ( Han et all l2013dY Our 
proof does not invoke the a priori upper bound of the departure rates (path flows). 

Keywords: dynamic user equilibrium; elastic demand; variational inequalities; differen- 
tial variational inequalities; existence 

1 Introduction 

This paper is concerned with an elastic de mand extens i on of t he fixed-demand dynamic traffic 
assignment model orig i nally presented in Friesz et al. (119931) and discussed subseq uently by 



Friesz and Mookherieel (|200d ); lFriesz et al.1 (|201 ll . 120131 ) ; iFriesz and Meimandl (|2013l ). As such, 
it is concerned with a specific type of dynamic traffic assignment known as dynamic user 
equilibrium (DUE) for which travel cost, including delay as well as early and late arrival 
penalties, are equilibrated and demand is determined endogenously. 



1.1 Introductory remarks on dynamic user equilibrium 

In the past two decades there have been many efforts to develop a theoretically sound for- 
mulation of dynam ic user equilibrium especially in continuous time. As is pointed out by 
Friesz et al.1 (j201ll h DUE models tend to be comprised of five essential components: 



1. a model of path delay; 
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2. flow dynamics; 



3. flow propagation constraints; 



4. a route and departure time choice model; and 

5. a model of demand growth (evolution). 

Dynamic user equilibrium models from the early 1990s forward have been largely concerned 
with the sub-models 1 through 4 above, concentrating on the within-day time scale for which 
travelers make routing and departure time decisions. Items 1 through 3 above form a sub- 
model known as the dynamic network loading (DNL), which determines arc- and/or path- 
specific volumes and flows as well as experienced path delays when departure rates (path 
flows) are known for each path. Item 4 aims at expressing mathematically the notion of Nash- 
like equilibrium conditions. Item 5, demand evolution, occurs on the day-to-day time scale 
and allows travel demands to be updated. 

There are multiple means of expressing the Nash-like notion of a dynamic user equilibrium, 
including the following: 



(i) an in finite-dimensional variational inequality (jFriesz et al.l . ll993l ; ISmith and Wistenl . ll994L 
19951 ) 



(ii) a nonlinear complementarity problem ( Wie et al. . 20021 ; Han et al. . 2011 ) 



(iii) a diff erential variational inequality ( Friesz et al. . 2001 . 20131 . 2011 ; Friesz and Mookheriee 
200d ) 



(iv) a differential complementarity system (jPang et all 12011 



The variational inequality representation is presently the primary mathematical form em- 
ployed for dynamic user equilibrium. One of the advantages of such variational inequality 
formulation is that it subsumes almost all simultaneous route-and-departure time choice DUE 
models regardless of the dynamic network loading procedure employed. In fact, the arc dy- 
namics and flow propagation constraints can be naturally embedded in the effective delay 
operator which is viewed as a mapping between two Hilbert spaces. However, it would be a 
mistake to think that somehow the VI formulation is an "easier" formulation since the effective 
delay operators are generally not knowable in closed form; in fact, the delay operators may be 
non-analytic and may need to be derived from an embedded delay model, data combined with 
response surface methodology, or data combined with inverse modeling. Analytical results of 
the effective delay operator is crucial for understanding qualitative properties of dynamic user 
equilibrium such as existence, uniqueness and convergence of certain computational schemes. 
The existence of DUE requires that the delay oper ator is continuo us, which is a consequence of 
a generalization of Brouwer's fixed point theorem ( Browder . 1968). In addition, the uniqueness 
of a D UE solution further req u ires t hat the delay operator is strongly monotone (iNagurneyi . 
19931 ). iFriesz and Mookheried l|2006l ) propose and test a fixed point algorithm implemented 
i n continuou s time to solve the differential variational inequality (DVI) formulation of DUE 
dFriesz et all l200lf >: ;hat algorithm requires monotonic effective delay operators to assure 
convergence. 
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1.2 Some history of DUE with elastic demand 



In this section we review some of the f ew prior efforts to model DUE with elastic travel 
demand. This review is based in part on Friesz and Meimand ( 20131 ). Most of the studies of 
DUE reported in the DTA literature are about dynamic user equilibrium with constant travel 
demand for each origin-destination pair. We denote such a dynamic user equilibrium with 
fixed demand by F-DUE. It is, of course, not gener ally true that travel dem and is fixed, even 
for short time horizons. Arnott et al. ( 19931 ) and Yang and Huang ( 19971 ) directly consider 
elastic travel demand; their work possesses a limited relationship to the analysis presented 
in this paper, for their work is concerned with a simple bottleneck instead of a nontrivial 
net work, which is our focu s . 

Yang and Meng (1 19981 ) extend a simple bottleneck model to a general queuing network 



with known elastic demand functions for each origin-destination (OD) pair. IWie et al.l (|2002l ) 
study a version of dynamic user equilibrium with elastic demand, usin g a complementarity 
formulation that requires path delays to be expressible in closed form. ISzeto and Lol ( 2004 ) 
study dynamic user equilibrium with elastic travel demand when network loading is based 
on the cell transmission model (CTM); their formulation is based on discr ete time and is 
expressed as a finite-dimensional variational inequality (VI) . » 

study dynamic 

user equilibrium with elastic travel demand for a network whose traffic flows are also described 
by CTM. 

Although liriesz et al.1 (|201ll ) show that analysis and computation of dynamic user equilib- 
rium with constant travel demand is tremendously simplified by stating it as a a differential 
variational inequality (D VI), they do not discuss how elastic demand may be accommodated 
within a DVI framework. Friesz and Meimand ( 20131 ) extend the DVI formulation to an elas- 
tic demand setting. Such a DVI formulation of elastic demand DUE is not straightforward. 
In particular, the DVI presented therein has both infinite-dimensional and finite-dimensional 
terms; moreover, for any given origin-destination pair, inverse travel demand corresponding 
to a dynamic user equilibrium depends on the terminal value of a state variable representing 
cumulative departures. The DVI formulation achieved in that paper is significant because it 
allows the still emerging theory of differential variational inequalities to be employed for the 
analysis and computation of solutions of the elastic-demand DUE problem when simultaneous 
departure time and route choice are within the purview of users, all of which constitutes a 
foundation problem within the field of dynamic traffic assignment. 

A good review of recent insights into abstract differential variational inequality theory, 
i nclud ing computational methods for solving such problems, is provided bv lPang and Stewart 
( 20081 ) . Also, differential variational inequalities i nvolving the kind of explicit, agent-specific 
control variables employed herein are presented in iFrieszl poiol ). 



1.3 Main results 



This 



(201 



aper further advances the knowledge of continuous-time E-DUE based on lFriesz and Meimand 
in terms of formulation and qualitative properties. In particular, we state and prove 



an equivalent variational inequality (VI) formulation of the E-DUE problem. The proof uses 
a measure-theoretic argument, and do es not invoke the opt i mal c ontrol framework which is 
the primary methodology employed by iFriesz and Meimand (l2013h to establish the DVI for- 
mulation. Such a result is presented in Theorem 13.11 

As an immediate application of the proposed variational inequality formulati on, we will 



analyze and establish existence result for the E-DUE problem. As commented in lHan et al.1 
(j2013ch . ;he most obvious approach to establishing existence is to convert the problem to an 
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equivalent variational inequality problem and then apply a version of Brouwer's fixed point 
existence theorem. Nearly all proofs of DUE existence employ Brouwer's fixed point theorem, 
ei ther impl i citly o r explicitly. One statement of Brouwer's theorem appears as Theorem 2 
of Browder ( 19681 ). Approaches based on Brouwer's theorem require the set of feasible path 



flows under consideration to be compact and convex in a topological vector space, and typically 
inv olve the a pr i ori bo u nd on path flows. For inst ance, using the link delay model introduced 
by iFriesz et al.1 (jl993h . Izhu and Marcottel (200(3) show that a route choice (RC) dynamic 



user equilibrium exists under certain regularity conditions. In their modeling framework, the 
departure rate at each origin is given as a priori and assumed to be bounded from abov e . Thu s 
one is assured that all path flows are automatically uniformly bounded. In Wie et al. ( 20021 ). 



the existence of an arc-based user equilibrium is established under the assumption that the 
path flows are a priori bounded. 

Difficulties arise in the proof of a general existence theorem from two aspects: (i) in a 
continuous-time setting, the set of feasible path flows is not compact; and (ii) the assumption 
of a priori boundedness of path flows, which is usually required for a topological argument, 
does not arise from any behavioral argument or theory. The existence proof provided by 
this paper manages to overcome these two major difficulties. Regarding item (i) above, we 
employ successive finite-dimensional approximations of the feasible path flows set, which allows 
Brouwer's fixed point theorem to be applied. Regarding item (ii), we propose an in-depth 
analysis and computation of the path flows under minor assumptions on the travelers' disutility 
functions. 

Existence result for the elastic demand case is further complicated by the fact that the 
corresponding VI formulation has both infinite-dimensional an d finite-d i mensi onal terms (see 
Theorem 13.11 below). In order to apply Browder's theorem ( Browder . 19681 ). one needs to 



work in an extended Hilbert space that is a product of an infinite-dimensional space and a 
finite-dimensional space, and define appropriate inner product that allows compactness and 
weak topology to be properly defined. It is significant that our existence result for E-DUE, 
stated and proved in Theorem 15.61 does not rely on the a priori upper bound of path flows 
and can be established for any dynamic network loading sub-model with reasonable and weak 
regularity conditions. 

1.4 Organization 

The rest of this paper is organized as follows. Section [2] recaps several key notations and 
concepts that are essential for the subsequent discussion. In Section O we present one of 
the main results of this paper: the variational inequality (VI) formulation of the simulta- 
neous route-and-departure choice dynamic user equilibrium with elastic travel demand (E- 
DUE) . Section [H briefly rey iews the differential variational inequality formulation proposed by 
Friesz and Meimandl ( 20131 ). Section [5] establishes the existence result for E-DUE based on 



the VI formulation proposed in Section The proof overcomes several difficulties known to 
other researchers, including non-compactness of the feasible set and the a priori boundedness 
of path flows. 

2 Notation and essential background 

The time interval of analysis is a single commuting period or "day" expressed as 

[*o,*/] cS 1 
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where tf > to, and both to and tf are fixed. Here, as in all DUE modeling, the single most 
crucial ingredient is the path delay operator, which provides the delay on any path p per unit 
of flow departing from the origin of that path; it is denoted by 

D p (t, h) Vp(£V 

where V is the set of all paths employed by travelers, t denotes departure time, and h is a 
vector of departure rates. From these, we construct effective unit path delay operators ^ p (t, h) 
by adding the so-called schedule delay / (t + D p (t, h) — Ta); that is 

%(t,h)=D p (t,h) + f(t + D p (t,h)-T A ) VptV (2.1) 

where Ta is the desired arrival time and Ta < tf. The function /(•) assesses a penalty 
whenever 

t + D p (t, h) £ T A (2.2) 

since t + D p (t, h) is the clock time at which departing traffic arrives at the destination of path 
p G V. We stipulate that each 

%(-,h) : [to,*/] — ► ^V^V 

is measurable and strictly positive. We employ the obvious notation 

[%(-,h):peV) g5? |p| 

to express the complete vector of effective delay operators. 

It is now well known that path delay operators may be obtained from an embedded delay 
model, data combined with response surface methodology, or data combined with inverse 
modeling. Unfortunately, regardless of how derived, realistic path delay operators do not 
possess the desirable property of monotonicity; they may also be non-differentiable. 



2.1 Review on dynamic user equilibrium with fixed demand (F-DUE) 

For the completeness of our presentation, in t his section we reca p the definition of DUE with 
fixed travel demand, originally articulated bv lFriesz et al. I (|l993h . 

Let us introduce the fixed trip matrix (Qij : (i, j) G W), where each Qtj G is the 
fixed travel demand between origin-destination pair G W, where W is the set of origin- 
destination pairs. Note that Qij represents traffic volume, not flow. Finally we let Vij C V be 
the set of paths connecting origin-destination pair € W. As mentioned earlier, h is the 
vector of path flows h = {h p : p G V}. We denote the space of square integrable functions on 
the real interval [to, */] by L 2 [to, tf]. We stipulate that each path flow is square integrable, 
that is 

he{L\[t Q , t f ]f l 

where (L+[to, */]) is the positive cone of the |'P[-fold product of the Hilbert space L 2 [to, tf]. 

Each element h = (h p : p G V) G (L^_[to, tf]) is interpreted as a vector of departure-time 
densities, or more simply path flows, measured at the entrance of the first arc of the relevant 
path. It will be seen that these departure time densities are defined only up to a set of measure 
zero. With this in mind, let v denote a Lebesgue measure on [to, tf], and for each measurable 
set S C [t , tf], let Vjy(t G S) denote the phrase for v-almost all t G S. If S = [to, tf], then we 
may at times simply write Vj,(t). 
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We write the flow conservation constraints as 



I h p (t) dt = Qij V(i, j) G W (2.3) 



to 



where (j2.3|) consists of Lebesgue integrals. Using the notation and concepts we have men- 
tioned, the feasible region for DUE when effective delay operators are known is 



An 



\h>0: [ f h p (t)dt = Q lJ V(i,j) G w\ C (L 2 + [t ,t f ]) lVl (2.4) 



In order to define an appropriate concept of minimum travel costs in the present context, we 
require the measure-theoretic analog of the infimum of a set of numbers. In particular, for 
any measurable set S C [i , tf] with v(S) > 0, and any measurable function / : S — >• 3?, the 
essential infimum of /(•) on S is given by 

essinf{/(s) : s e S} = sup{xG$R: u{s £ S : /(a) < x} = 0} (2.5) 

Note that for each x > essinf{/(s) : s € S} it must be true by definition that 

v{s£S : f(s) < x} > 

Let us define the essential infimum of effective travel delays 

v p = v p (h) = essinf{* p (t, h) : t G [in, tf]} > Vp G V 

Vij = Vij(h) = min{v p (h) : p G Vij} V(i,j)£W 



The f ollowing definition of dynamic user equilibrium was first articulated by iFriesz et al- 
ii 19931 ): 

Definition 2.1. (Dynamic user equilibrium) A vector of departure rates (path flows) 
h* G An is a dynamic user equilibrium if 

h* p (t) > 0, p G Vij => %(t, h*) = VlJ (h*) y u (t) G [t , t f ], V(i, j) G W (2.6) 

* p (t, h*) > Vij(h*) y u (t) G [t , t f ], V(i, j) G W (2.7) 

H^e denote this equilibrium by DUE{^, Ao, [toji/])- 

Using measure theoretic arguments, Friesz et al.l (|l993l ) establish that a dynamic user 
equilibrium is equivalent to the following variational inequality under suitable regularity con- 
ditions: 

find h* G Ao such that 
E [ tf * P (t,h*)(hp-h*p)dt>0 }VI(V,A ,[to,tf]) (2i 
Vh G A 



1 The demand satisfaction expressed via Lebesgue integrals is not enough to assure that the path flows are 
bounded. This has been the major hurdle to proving existence without the a priori invocation of bounds on 
path flows. 
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2.2 Dynamic user equilibrium with elastic demand (E-DUE) 

The general setup of DUE with elastic demand is similar to that of the fixed demand case, with 
the exception that the total travel demand Qij between an origin-destination pair (i, j) G W 
is no longer a prescribed constant. Rather, transportation demand is assumed to be expressed 
as the following invertible function 

Qij = Fij[v] 

for each origin-destination pair (i, j) G W, where W is the set of all origin-destination pairs and 
v is a concatenation of origin-destination minimum travel cost associated with (i, j) G W. 
That is, we have that 

G 

v = (vij : (i, j) € W) G »f 1 

Note that to say is a minimum travel cost means it is the minimum cost for all departure 
time choices and all route choices pertinent to origin-destination pair G W. Further note 
that Qij is the unknown cumulative travel demand between G W that must ultimately 
arrive by time tf. 

We will also find it convenient to form the complete vector of travel demands by concate- 
nating the origin-specific travel demands to obtain 

Q = (Qij ■ (i, j) G W) G 5?f 
F: ^l^sftfl, „ Q 
The inverse demand function for every G W is 

%• = &ij [Q] 

and we naturally define 

Q[Q] = {€> ij [Q]:(i,j)€W)€3lW 

As a consequence, we employ the following set of feasible departure flows when the travel 
demand between each origin-destination pair is unknown. 



A 



(h, Q) : h>0, f f h p (t)dt = Qij V(i, j) G W \ C {L 2 [t , t f ]) |P| x ^ l 



(2-9) 

where (-L 2 [^o 5 Vl)' 7 '' x ^4^' ^ t ne direct product of the l^-fold product of Hilbert spaces 
consisting of square-integrable path flows, and the |W|-dimensional Euclidean space consisting 
of vectors of elastic travel demands. 

With preceding preparation, we are in a place where the simultaneous route-and-departure- 
choice dynamic user equilibrium with elastic demand can be rigorously defined, as follows. 

Definition 2.2. (Dynamic user equilibrium with elastic demand) A pair (h* , Q*) G A ; 
where h* is a vector of departure rates (path flows) and Q* is the associated vector of travel 
demands, is said to be a dynamic user equilibrium with elastic demand if for all (i, j) G W, 

h* p (t) > 0, P eVij => %(t,h*) = @ij[Q*} V„(t) g [to, t f ] (2.10) 

%(t, h*) > &ij[Q*} V„(t) G [t , */], VP e Vij (2.11) 
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3 The variational inequality formulation of E-DUE 

Experience with differential games suggests that the DUE problem with elastic demand can 
be expressed as a variational inequality, as shown in the theorem below. 

Theorem 3.1. (E-DUE equivalent to a variational inequality) Assume VP p (-, h) : 
[to, tf] —> is measurable and strictly positive for allp £? and all h such that (h, Q) G A. 
Also assume that the elastic travel demand function is invertible with inverse Oy[-] for all 
(i, j) G W. Then a pair, {h* , Q*) G A, is a DUE with elastic demand as in Definition \2.£\ if 
and only it solves the following variational inequality: 

find (h* , Q*) G A such that "1 

£ [ f %(t,h*)(h p -h;)dt- Y 9ij[Q*][Q»-Qij\ > \vi(%e,to,t f ) (3.12) 

V(h, Q) G A 

Proof, (i) [ Necessity] Given a DUE solution with elastic demand (h* , Q*) G A, we easily deduce 
from ([27HJ|) and ^JT\i that for any (h, Q) G A, 

Y h*){h P - h* p )dt - Y, % [Q*\ [Qij - Qij] 

= E ( E f f %(t,hnh p (t)dt-e ij [Q*]-Q ij \ 

- E E f *f (*> dt - e ^ ■ Qh 

= e ( e r ^)^(*) * - 9i i^i • ^ j 

- E ( E e * • W dt - e « [^1 • Q*i 



E ( E /*' *p(*> ^M*) * - e ^*l • G« ) 

- e wi- ( e r 

(i,i)eW \pe-Pij *° 



/i*(t) - Q| 
'to 



E ( E /*' *p(*> * - e *i[Q*] • o« ) ( 3 - 13 ) 
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Observe that in (|3.13p . given any (i, j) G W, 

%(t, h*)h p (t) dt - OijiQ*] ■ Qij 



E 



to 



> v*j Y !' hp(t) dt - QijlQ*} ■ 
peVn Jt ° 




h p (t) dt - Qi 



\ 



(3.14) 



/ 



where v*j is the essential infimum of ^ p (-, h*) on [to, tf] for all p G Vij, and is equal to @ij[Q*] 
according to (]2.10p . 

As an immediate consequence of (|3. 13[) and (|3.14p . the following inequality holds for all 
(h, Q) G A. 

Y I 1 9 p (t, h*){h p - h* p )dt - Y % IQ*) [QH - Q*j] > ( 3 - 15 ) 



pev 



(i,j)€W 



(ii)[Sufficiency] Assume that (]3. 15|) holds for any (h, Q) G A, then h* is a solution to the 
DUE problem with fixed demand given by Q*, since in this case the second term in the left 
hand side of f|3. 15j) vanishes and we recover the well-known VI for the fixed demand case, see 
([Z5]) . Therefore by Definition EH for any (i, j) G W, 



V„(t) G [t , t f ] 



| h* (t) > 0, p G Vij * p [t, h* (t)] = «* 

\* p (t, fc*) > 4 ft*(t) = o 

In order to show that (h*, Q*) is a DUE with elastic demand using definition (|2.10p and (|2,lip . 
it suffices to establish that v*- = 0y[Q*], V(z, j) G W. We proceed as follows. 
Fix arbitrary (k, I) G W such that Q* kl ^0@. We define the pair (h, Q) G A as 



ah*(t) 

h;(t) 

a Qh 
Q*j 



p£V k i 
peV\V k i 

(i, J) = (*, 

(i,i)€W\(fc, 



Vt G [t , t f ] 



where a G is an arbitrary positive parameter. Substituting (h, Q) for (h, Q), the left 
hand side of (|3.15p becomes 

rtf 



p£V 



to 



* p (t,h*)(h p -h* p )dt- Y ®ijiQ*] 

(i,i)ew 



^ [ f y p (t,h*){ah;-h* p )dt+ Y 



h*){hl-ht)dt 



®ki[Q*} [aQh-Ql 



(i,j)eW\(fc,0 



(a - 1) Y ^ ^) W * - (a - l)©fci[Q*] Qfci 



(a-lK^-(a-l)e w [Q*]Q^ 



2 For any origin-destination pair (i, j) such that Q*j — 0, the logical statements (|2,10|l and (|2.11[) automati- 
cally hold, since in this case hp(-) vanishes almost everywhere for all p £ Vij. 
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We conclude from (|3.15|) that 



(a-l)K,-e w [Q*])C& > 

Since a G is arbitrary, there must hold t;^ = O^fQ*], for any (k, I) G W. The proof is 
complete. □ 



4 The differential variational inequality formulation of E-DUE 

The notion of dynamic user equilibrium can be alternatively illustrat ed and analyzed using th e 
mathematical paradigm of differential variational inequality (DVI) ( Friesz et al. . 2001 . 201 ll ). 
In particular, the demand satisfaction can be easily rewritten as a two-point boundary value 
problem 



dy 



ij 



dt 



h p&> 



Vij(tf) = Qij V(i, j) G W 



(4.16) 



where yij(-) is interpreted as the cumulative departure curve. In the context of the elastic 
demand case, we assume there are unknown terminal states yij{tf), for all 6 W, which 
are the realized DUE travel demands. Moreover, for each origin-destination pair € W, 
inverse travel demand is expressed as 



%■ = & ij [y (if)} 



(4.17) 



where 



y = ( yij (-): (i,j)EW)e(H 1 ([t ,t f })) 



\w\ 



'H 1 {[tQ, tf]) is the Sobolev space consisting of weakly differentiable functions whose weak 
derivatives are square- integr able. Thus, the yij(tf), for all G W, will be determined 

endogenously to the differential variational inequality pres ented subseq u ently in Theorem 14.11 
Such an approach contrasts to the approach employed by iFriesz et al.l ()201lh to study fixed- 
demand DUE by making each yij(tf) an a priori fixed constant Qij. Accordingly, we introduce 
the following dynamics: 



dyij 
dt 



yy(to) = V(ij)eW 



(4.18) 



As a consequence, we employ the following alternative form of the feasible set: 



Ai = < h > : 



dt 



Vijih) = V(i,j) G W \ C (L\ [t ,t f ]) lVl 



pePi, 



(4.19) 



Note that the feasible set K\ in ()4.19j) is expressed as a set of path flows since knowl- 
edge of h completely determines the demands that satisfy the initial value problem (|4.16p . 
The differential vari a tiona l inequality formulation of E-DUE is first stated and proved by 
Friesz and Meimandl ( 20131 ) using optimal control theory. 



Theorem 4.1. (E-DUE equivalent to a differential variational inequality) Assume 
*&p(-,h) : [to,tt] — > ^++ * s measurable and strictly positive for all p G V and all h G Ai. 
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Also assume that the elastic travel demand function is invertible, with inverse (Q) for all 
G W. A vector of departure rates (path flows) h* £ Aj is a dynamic user equilibrium 
with associated demand Q*(tf) if and only if h* solves the following DVI(*f/,@,to,tf). 

find h* G Ai such that } 
[ 1 %(t, h*)(h p - h* p )dt - ^ [y* (t f )] [ yij (t f ) - y*j (t f )] > I W J(*, , t , t f ) 



(ij)6W 

V/i g Ai 



Proof. See Friesz and Meimand ( 20131 ). 



(4.20) 
□ 



5 Existence of dynamic user equilibrium with elastic demand 

In this section, we will establish existence result for VI(^>, 0, to,tf), an equivalent formulation 
of E-DUE. Our proposed approach is meant to incorporate the most general dynamic network 
loading sub-model with minimum regularity requirements, and to yield existence of E-DUE 
without invoking the a priori upper bound on path flows. In subsequent analysis, we will 
rewrite VI(^, 0, to, tf) as a variational inequality problem in an extend e d Hil bert space and 
then employ a version of the Brouwer's Fixed Point Theorem (Browder, 19681 ). To this end, 
we introduce the product space E= (L 2 [t , t f ]) lPl x RlWI, which is a Hilbert space with the 
induced inner product defined as follows 

\-p\ r t f \w\ 
(X,Y) E ± Y, m-m(t)dt + Y^j (5-21) 

i=l Jt ° j=l 

X = (6(0, ...,^ v \{-),u 1 , ...,u M ) G E (5.22) 
Y = (m(-), • • • , V\v\(-), v 1} . . . , v m ) G E (5.23) 

Let us recall the VI formulation presented in Section EJ the set A of admissible pair (h, Q) 
can now be embedded in the extended space E: 

A= (h, Q) G (Lllt , t f ])M x »™ : f f h P (t)dt = Q tJ V(i,i)eW CE 

pep* Jt ° J 

In view of the inverse demand function B = (Gy : (i, j) G W), we introduce notation 

6- = (-Qij : (i, j) G W) : ^f 1 — »• kL W| 
We next define the mapping 

T : A — > E, (h,Q) ^ (*(•, h), e~(Q)) (5.24) 

where (h, Q) G A, *(-, h) G (L* [t , */]) |W| , G »!^'- Such a mapping is clearly well- 

defined. With the preceding discussion, the VI formulation of the DUE problem with elastic 
demand is readily rewritten as the following infinite-dimensional variational inequality in the 
extended Hilbert space. 
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find X* £ A such that 
(T(X*), X -X*) E > 
VlGA 



} VI(T,e,A 



to,*/) 



(5.25) 



where X = (h, Q) and X* = (h* , Q*). Problem (|5.25|) is written in the generic form of 
variational inequality, wh ich allow s analy sis regarding solution existence to be carried out in 
a framework provided by Browder ( 19681 ). 



5.1 Existence result for VI (F, 6, A, t , tA 



Our qualitative analysis regarding solution existence for the variational inequality (|5.25[) is 
based on the following extension of Brouwer's fixed point theorem to topological vector spaces. 



Theorem 5.1. (Browdei . 1964 ) Let K be a compact convex subset of the locally convex 
topological vector space V, T a continuous (single-valued) mapping of K into V* , where V* 
is the dual space of V . Then there exits uq in K such that 



(T(uo), uq-u) > 



for all u £ K . 



Proof. See Browder ( 19681 ) . 



□ 



In preparation for our existence proof, we recap several key results from functional analysis 
that facilit ate our presentation. In part icular, we note the following facts. The reader is 
referred to Royden and Fitzpatrick ( 19881 ) for more detailed discussion on these subjects. 



Proposition 5.2. The space of square-integrable real-valued functions on a compact interval 
[to, tf], denoted by L 2 [to, tj], is a locally convex topological vector space. In addition, the 

is also a locally convex topological 



\v\ 



\V\-fold product of these spaces, denoted by {l?\t§, tff\ 
vector space. 

Proposition 5.3. The dual space o/(L 2 [io, t/])'^' has a natural isomorphism with (L 2 [io, i/])'^'. 
The dual space of the Euclidean space consisting of columns of |W| real numbers is inter- 
preted as the space consisting of rows of |W| real numbers. As a consequence, the dual space 



w 



of (L 2 [t , tf}) 1 ' 1 x 5?l w l is again (L 2 [t , t/]) 1 ' 1 x 5? 

Proposition 5.4. In a metric space (therefore topological vector space), the notion of com- 
pactness is equivalent to the notion of sequential compactness, that is, every infinite sequence 
has a convergent subsequence. 

Theorem l5.1l is immediately applicable for showing that Vl(F, 0, A, to, tA has a solution 
if (1) T is continuous; and (2) A C E is compact. Unfortunately, such compactness does not 
generally hold for the problem we study herein. To overcome such an obstacle, we proceed 
in a similar way as in lHan et al.l ()2013d ) by considering finite-dimensional approximations 
of the underlying infinite-dimensional Hilbert space. Another major hurdle that stymied 
many researchers is the a priori upper bound of path flows. Such bound is important for a 
topological argument that we will rely on in the proof, but does not arise from any physical 
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or behavioral perspective of traffic modeling. In fact, as observed by Bressan and Han ( 201ll ). 
the equilibrium path flows could very well become unbounded or even produce dirac-delta, if 
no additional assumptions are made regarding exogenous parameters of the Nash-like game, 
such as travelers' disutility functions. 

The following assumptions are key to our analysis of path flows. The first assump- 
tion, (Al), poses hypothesis on drivers' perceived arrival costs; the second assumption, 
(A2), is concerned with the model of link /path dynamic s and can be easi ly satisfied by 
existing mo dels such as the V ickrey model ( Vickrev . 19691 ; Han et al. . 2013a b ), the LWR- 
Lax model (IFriesz et all 12013!) . and the Lighthill-Whitham- Richards model (IDaganzd . 11994 : 
Lighthill and Whithaml . Il95a iRichardl Il95d ) . 



(Al). The function /(•) appearing in ()2.ip is continuous on [to, tf] and satisfies 
/(* 2 )-/(ii) > A(i 2 -*i) Vt < h < t 2 < t f 

for some A > — 1 



(5.26) 



(A2). The first-in- first- out (FIFO) rule is obeyed on a path level. In addition, each link 
a € A in the network has a finite exit flow capacity M a < oo. 



Remark 5.5. Assumption (Al) is employed in an in- depth analy s is of th e network flow later 
in the proof of Theorem 15.61 The reader is referred to Han et al. 1 2013a) for the motivation 
and generality of such an assumption. Notice that if /(•) is continuously differentiable, then 
(Al) is equivalent to requiring that f'(t) > —1, t € [to, tf]. 

With the preceding preparation, we are now ready to state and prove the main result of this 
section. 

Theorem 5.6. Assume that the effective delay operator ^ : Ao — > (L 2 [to, t/]) is continuous. 
In addition, let assumptions (Al) and (A2) hold. If the inverse demand function G : — > 
is continuous, then the variational inequality Vl(F, G, A, to, tA has a solution. 



Proof. Given that \& and G are both continuous mappings, it is straightforward to verify by 
definition (|5.24j) that the mapping T : A — > E is also continuous. 

Since Theorem 15 . 1 1 cannot be directly ap plied to obtain a so lution of the infinite-dimensional 
VI, let us instead employ the technique from lHan et al. ( 2013cJ ) by considering finite-dimensional 
approximations of A. More specifically, consider for each n > 1 the uniform partition of [to, tf] 
by n sub- intervals J , . . . , I n . Define the finite-dimensional subset of A: 



Ar, 



(7n (•),..., VlO> Qu---,Q\w\) eA: 



hi(-) is constant on each P VI < j < n, VI < i < \V\ > (5.27) 



Moreover, it is not restrictive to assume that there is an upper bound on the elastic demand 
for each origin-destination pair. That is, there exists a vector U = (Uij) € 3?!^ such that 



< Qij < V(i,j)eW, V(/i, Q)gA 
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We can then show that A n defined as such is convex and compact in A for each n > 1. The 
proof is postponed until Lemma 17.11 of the Appendix for the clarity and concision of our 
presentation. 

We are now in a position where Theorem 15 . 1 1 applies to each A n . In other words, fix n > 1, 
there exists some X n '* = (h n '*, Q n '*) G A„ such that 

{T{X n '*), X n - X n '*) E >0 VI n G A„ (5.28) 

As a consequence of (|5.28|) . we have 

52 f f %{t,h n >*)h n '*(t)dt- 52 & iAQ n nQ r :f 



< 



52 h n >*)h n (t)dt- 52 ®iAQ n '*m ( 5 - 29 ) 



for all (h n , Q n ) G A„. In particular, for all h n G (Z/j_[i , i/])'^' such that 

rtf 



inequality (|5.29|) becomes 



52 I* K(t)dt = Q«f v(i,i)ew 



52 I' *p(*> h n '*)h n '*(t)dt < 52 f* %(t, h n '*)h n (t)dt (5.30) 

Recall that h n '*(-) is piecewise constant, thus (I5.30P implies the following: 

C*(t) > 0, iG/ fc => / %(t,h n '*)dt = min / ^Jt,h n '*)dt (5.31) 

for all p G (i, j) G W. We invoke Lemma 17.21 from the Appendix and find a constant 
M G K ++ such that 

/£'*(*) < M VtG[t , */], VpGP, Vn>l 

With the uniform upper bound .M on the path flows h n '* and upper bound Uij for each Qjj, 
(i, j) G W, one can extract a subsequence {X nfe '*} out of {X n '*} such that X nk '* converges 
weakly to some X* G A as k — > oo, where the weak topology on £" is defined via the inner 
product (■, ■) E . 

In view of (|5.28p . by strong continuity of J- and the weak convergence X n> * —> X* , we 
conclude that 

(T(X*), X — X*) E > VXGA 

□ 

Remark 5.7. The existence of an E-DUE requires that the embedded network loading pro- 
cedure yields a delay operator tha t is contin u ous. S uch regularity condition coincides with 
that of the fixed demand case, see \Han et al. Theorem 15.61 subsumes all notions of 

simultaneous-route- and- departure choice user equilibrium regardless of the arc dynamic, flow 
propagation and delay model employed. 
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6 Conclusion 



We have shown that dynamic network user equilibrium based on simultaneous departure time 
and route choice in the presence of elastic travel demand may be formulated as a variational 
inequality (VI)- Su ch a result is a nontrivial extension of the VI formulation established by 
Friesz et all (| 19931 ) for dynamic user equilibrium with fixed travel demand. The proof relies 



on a measure-theoretic argument and provides unique insights on the qualitative properties of 
E-DUE such as existence, which is also established in this paper. 

Similar to the fixed demand case, existence of the elastic demand case is easily analyzed 
using the framework o f Brouwer's fix ed point theorem for infinite-dimensional variational 
inequalities proposed in Browder ( 19681 ). Nevertheless, the variational inequality for the elastic 
demand case has both infinite- and finite-dimensional terms. The proof of existence requires 
construction of an extended space that subsumes both parts, as well as an appropriate choice 
of inner product that allows compactness and weak topology to be defined. It is significant 
that our proof does not rely on the a priori upper bound on path flows. 

We are conducting further investigation of the E-DUE problem in terms of computation, 
based on the DVI formulation. In particular, the optimal control framework inherent in the 
DVI formulation allows us to device an iterative scheme based on a fixed-point-problem refor- 
mulation. Convergence of such a scheme requires some sort o f monotonicity to be articulated 
and proved, such as those mentioned in iFriesz et al.1 (Il993h . which is the focus of ongoing 
research. 
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7 Appendix 

Lemma 7.1. Assume that there exists a vector U = (£%) G 3?|^ such that 

< Qij < Uij V(i, j) G W, V(/i, Q)gA 
Then for each n > 1, the subset 

A n = j(/ii(-),..., h\ v \{-), Qi,...,Q m ) G A: 

hi(-) is constant on each P VI < j < n, VI < i < \V\^j (7.32) 

is compact and convex in A. 

Proof. We begin with verifying convexity. Let X = (h, Q) and X = (h, Q) be any two 
elements of A n . Given any a € (0, 1), we have that 

I 1 (ah p (t) dt + (1 - a)h p (t)) dt = a Q - + (1 - a) C? - V(i, j) G W 

Moreover, a h p (-) + (1 — a) h p (-) clearly remains constant on each sub-interval P , j = 1, . . . , n, 
for all p G V. We thus conclude that aX + (1 — a)X G A n . 

Next, let us investigate compactness. From now on let us fix n > 1. In light of Proposition 
15.41 it suffices to establish sequential compactness for A n . We consider an arbitrary infinite 
sequence {X fc } fc>1 C A n where X k = (h k , Q k ). For each k > 1 and p G V, let n k = (pi j) G 

be such that 

We then define p, k G Sff" to be the concatenation of all vectors fj, k , p G V . We also notice 
that the vectors p k , k > 1 are uniformly bounded by the constant 

n 

max Ui j ■ 

(»,i)ew t/-*o 
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Thus by the Bolzano- Weierstrass theorem, there exists a convergent subsequence {/• ffc '} fe /> 1 - 

It is immediately verifiable that the corresponding subsequence h k ' converges uniformly on 
[to, tf] and also in the L 2 norm. 

Moreover, by virtue of the uniform bounds {/«, (i, j) € W, there exists a further subse- 
quence C {£/} such that Q k is a convergent subsequence according to the Bolzano- 
Weierstrass theorem. Thus, the subsequence {X fe "} fc//>1 converges with respect to the norm 
induced by inner product (|5.2ip . □ 



Lemma 7.2. Assume that (Al) and (A2) hold. Then there exists Ai > such that for all 
(h n <*, Q n '*) satisfying (fOSD . there hold 

h™>*{t) < M Vt€[to, */], Vp£V, Vn>l (7.33) 

Proof. In view of (A2), we are prompted to define the following constant 

M max = maxM a < + oo 

where A is the set of links of the network. Recalling the constant A from (Al), we choose 
constant A4 such that 

3M max 

M > -r — 

A + 1 

We claim that (|7.33[) holds for such M. Otherwise, if (|7.33p fails, there must exist some 
m > 1, q € V and 1 < j < m such that 

h™>*(t) = x > M teP 

Without losing generality, we assume that j > 1 and consider the interval P . By possibly 
modifying the value of the function *& q (-, h m '*) at one point without changing the measure- 
theoretic nature of the problem, we obtain t* £ P~ 1 such that 

%(f, h m >*) = sup * g {t, h m >*) 

We denote by r q (t, h m '*) the time of arrival at destination of driver who departs at time t 
along path q. According to the first-in-first-out (FIFO) principle, we deduce that Vt € P , 



(t-t j )X < J h^*(t)dt < M max (T q (t,h m >*)-T q (t*,h m >*)) 

where t 3 is the left boundary of the interval P . We then have the following estimation: 

^ q (t, h m '*) - -$ q (t*, h m '*) (7.34) 
= D q (t, h m >*) + f (r q (t, h m >*) - T A ) - D q {t\ h m '*) - f(r q (t*, h m >*) - T A ) 
> r q (t, h m >*) - r q (t*, h m '*) -(t- t*) + A(r 9 (t, h m >*) - r q (t*, h m >*)) 
= (A + l)(r q (t, h m n - r q (f, h m n) ~ (t ~ t*) 

>(A + l)^(t-^)-(i-n VteP (7.35) 
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Integrating (|7.35|) with respect to t over interval P shows the following: 



[ Vgit, h™'*) dt - (t j+1 - t^ q {t*, h m >*) 



> ^-Pf (A + 1)A _ / _ 

where t J , t 1+1 are respectively the left and right boundary of P . Since t* G P , we have that 



With this observation, (|7.36j) becomes 



(fj+l_^-)2 /(A + 1)A 



2 V M r ' 

> 



This implies 

f $> q (t,h m '*)dt > [ $>(t,h m '*)dt 

Jp Jli- 1 
which yields contradiction to (|5,3ip . This finishes the proof. □ 
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